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ABSTRACT 



Context. Polarization carries information about the magnetic fields in interstellar clouds. The observations of polarized dust emission are used 
^ to study the role of magnetic fields in the evolution of molecular clouds and the initial phases of star formation. Therefore, it is important to 
i—^ understand how different cloud regions contribute to the observed polarized signal. 

C/3 _ Aims. We study the grain alignment with realistic simulations, assuming the radiative torques to be the main mechanism. The aim is to study 
Ci , the efficiency of the grain alignment as a function of cloud position and to study the observable consequences of these spatial variations. 

Methods. Our results are based on the analysis of model clouds derived from MHD simulations of super- Alfvenic magnetized turbulent flows. 
, The continuum radiative transfer problem is solved with Monte Carlo methods to estimate the three-dimensional distribution of dust emission 
^ ■ and the radiation field strength. The anisotropy of the radiation field is taken into account explicitly. We also examine the effect of grain growth 
in cores both to the observed polarization and to the inferred magnetic field. 
■ Results. Using the assumptions of Cho & Lazarian, our findings are generally consistent with their results. However, the anisotropy factor 
00 is lower than their assumption of y = 0.7, and thus radiative torques are less efficient. Compared with our previous paper, P/I relations are 
steeper. Without grain growth, the magnetic field of the cores is poorly recovered above a few Ay. If grain size is doubled, the polarized dust 
emission can trace the magnetic field lines possibly up to Ay ~ 10 magnitudes. However, many of the prestellar cores may be too young for 
grain coagulation to play a major role. The inclusion of direction-dependent radiative torque efficiency weakens the alignment. Even with 
doubled grain size, we would not expect to probe the magnetic field beyond a few magnitudes in Ay . 
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1. Introduction the interstellar magnetic field. However, to explain the grain 

alignment needed for the polarization, this mechanism requires 
Magnetic fields are important in many astrophysical processes, rnuch stronger magnetic fields than have been observed. Purcell 
In star formation, the magnetic pressure may support an oth- (| 19791 1 introduced several processes that could cause the grains 
erwise gravitationally unstable cloud core against collapse, to become very fast rotators, and suggested that H2-ejection 
Magnetohydrodynamical models allow predictions of what might be a major cause of fast grain rotation. However, fur- 
magnetic fields might look like in interstellar clouds and pro- ther investigations of Ha-ejection identified several related pro- 
tostellar cores, but the models need to be tested against ob- cesses that make it inefficient in aligning dust grains, such as 
servations. Information about interstellar magnetic fields can grain wobbling (e.g., Jones & Spitzer 1967; Lazarian [T994I 
be derived with different techniques. Amongst the most useful Lazarian & Roberge lTWTI l. grain flipping (Lazarian & Draine 
are the polarization of starlight from background stars and the |1999al , and "nuclear relaxation" (Lazarian & Draine 1 1999b] l. 
polarized thermal dust emission at longer wavelengths, which These mechanisms and others were discussed in a review pa- 
arise from the alignment of the spin axis of the dust grains per by Lazarian (|2007l l. 

along the magnetic field. ^j^^ radiative torque mechanism has become a strong can- 

The main question with dust polarization is the alignment didate for the primary mechanism of grain alignment inside 

mechanism, or rather, what makes the dust grains spin up molecular clouds. The radiative torque mechanism involves 

in the first place? One of the earliest mechanisms proposed the transfer of momentum by collisions of photons onto the 

was the paramagnetic mechanism (Davis & Greenstein 1951 ), grain, causing a torque that rotates the grain around its axis, 

which is based on the direct interaction of rotating grains with It was first introduced by Dolginov d 19721 1 and Dolginov & 
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Mytrophanov d 19761 1. The efficiency of the radiative torques 
has been demonstrated using numerical simulations (Draine & 
Weingartner 19961 11997 ) and in a laboratory setup (Abbas et 
al.|2004). The predictions of the radiative torque mechanisms 
are roughly consistent with observations (e.g., Lazarian et al. 
[1997 i Hildebrand et al. '2000^. 

However, Draine & Weingartner ( I1996I I were incorrect to 
assume that the grains would just be spun up by radiative 
torques, and aligned by paramagnetic relaxation. In their sub- 
sequent paper, Draine & Weingartner ( 119971 ) confirmed that 
radiative torques would be able to align grains on their own, 
without relying on paramagnetic relaxation. Lazarian & Draine 
(I1999al l showed that thermal flipping would prevent the para- 
magnetic alignment mechanism being eff'ective in small grains. 

Cho & Lazarian (I2005I I used a spherically symmetric model 
cloud to calculate the size of grains aligned by radiative 
torques, assuming a constant anisotropy factor, y - 0.7, and 
neglecting the isotropic component of the radiation. Their cal- 
culations showed that even deep inside GMCs (Ay < 10), large 
grains can still be aligned by radiative torques. Cho & Lazarian 
(I2005i) presented an empirical formula for the minimum size 
of the aligned grain a^ig as a function of the density nu and 
extinction Ay, which is 



flaig = {log ne) (Ay.iD + 5)/2800;um. 



(1) 



The formula is, strictly speaking, only valid for spherically 
symmetric clouds, and carries the additional assumption of a 
constant anisotropy factor, y - 0.7. Pelkonen et al. (2007 ) used 
this formula to model polarization in a clumpy cloud. 

The obvious next step is to complete radiative transfer cal- 
culations in more detail, derive the actual anisotropy within an 
inhomogeneous cloud, and use it to calculate the efficiency of 
radiative torques. Similar work was carried out by Bethell et al. 
(I2007I I. These authors were able to derive f // relations remi- 
niscent of the observed ones, and their results agreed with the 
qualitative results of Cho & Lazarian (i2005i . They found that 
simplifying assumptions about the radiative anisotropy (e.g., a 
constant value of y) have only a moderate effect on the emer- 
gent degree of polarization. 

In the present paper, we study the anisotropy factor y and 
calculate the minimum sizes of the aligned grains within our 
model clouds. In addition to calculating ID models to compare 
with Eq. 1, we also compare our results with those of our pre- 
vious 3D models (Pelkonen et al. i2007l l. We investigate quanti- 
tatively the effects that grain growth has on the polarized signal 
expected from dense cloud cores. In particular, we examine a 
model where the grain growth is restricted to the densest re- 
gions. Furthermore, we examine where the observed polariza- 
tion actually originates along the line of sight. 

In the past few years, our understanding of the radiative 
torques has been improving rapidly. Lazarian & Hoang ( 20071 1 
presented the first analytic model of radiative torques, allowing 
for quantitative treatment of radiative torque alignment. At the 
time of submission of this paper, Hoang & Lazarian (.2008bl) 
expanded their treatment of radiative torques from beam ra- 
diation fields to dipole and quadrupole fields, and illustrated 
the dependence of radiative torque efficiency on the angle be- 
tween the magnetic field and the radiation. These important 
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Fig. 1. Polarization reduction factor as a function of Oaig 
using the normal {solid line) and doubled {dashed line) Li & 
Draine (1200 11 1 grain sizes. 
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Fig. 2. Radiative torque efficiency Qy for different grain sizes, 
fl, as a function of A (J. Cho, private communication). 



steps should be included in future studies of polarization within 
clumpy clouds. In this paper, we present an example where we 
take into account the angle between the local magnetic field 
and the direction of the anisotropic beam radiation to quantify 
the effect that it will have on emergent polarization. 



2. The basic equations 

To calculate the alignment efficiency of radiative torques, 
we follow the formalism presented in Draine & Weingartner 
(119961) . The rotational damping time due to gas drag is 



^drag.gas 



na^pa^^ .(8.74xl0Vr) 



35«H(27rmH/fcr)i/2 

ai SOOOcnr^ 

^—p^a^^T^ ( ), 

5 n^T 



(2) 



where ai is a geometric factor related to the moment of iner- 
tia (see Eq. 38 and Table 2 in Draine & Weingartner 1 19961 ). 
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Fig. 3. Av map of the high resolution 3D model as seen from the z-direction. The selected core regions are marked with white 
boxes. 



6 X ffi, p is the solid density, p3 = p/3 g cm"^, fleff is the ra- 
dius of a sphere of an equal volume, fl_5 = Oeff/lO"^ cm, ne is 
the number density of hydrogen atoms, T is the gas tempera- 
ture, and r2 = K. The rotational dampening time due to 
thermal emission of photons by a grain heated by starlight to a 
temperature is 



(1.60 X 10^ yr) 



(gabs) 18 K' U,,d 



(3) 



where misrf and M,ad are the energy density of the interstellar 
radiation field outside the cloud (Mathis et al. I1983I I and the 
radiation field illuminating the grain, respectively, and 



N 1 r 

ab.s; - M,) 

Mrad J 



QM dA, 



(4) 



and gabs ^ Ar ■ 

If we assume that the grain is illuminated by an isotropic 
radiation component m^™ and an anisotropic radiation compo- 
nent M™, a steady radiative torque gives the grain an angular 
velocity (Eq.65 in Draine & Weingartner .l996j 



xmGr>+<d<Qr)-ai)(- 



-), 



(5) 



where A is the average wavelength, oth the mass of a hydro- 
gen atom, (gp°) and (Qr) are the isotropic and anisotropic ra- 
diative torque efficiency, respectively, and ai is the unit vector 
of the rotational axis. However, this is an approximate treat- 
ment, since we ignore the dependence of radiative torque effi- 
ciency on the angle between the magnetic field and the radia- 
tion. Radiative torque efficiency is lower by more than a factor 
of ten for angles close to being perpendicular 

The randomization of the rotation of the grain is caused 
by collisions with gas molecules. When the grain rotates much 
more rapidly than the thermal rotation rate, 

' (6) 



(jJr - 



eff 



this randomization is greatly reduced. Thus, a suprathermally 
rotating grain is expected to be aligned with the magnetic 
field. The factor between the rotational rate due to the radiative 
torques, Wiad, and the thermal rotational rate, (jJt, is an indicator 
of the efficiency of the grain alignment, given by 

1 



^Wrad^2 _ _5ffl_paeff^ 

(jJt 1926^ niu nnkT 
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Fig. 4. The anisotropy factors at different wavelengths inside a 
ID model cloud. The thick lines are without and the thin lines 
with dust emission. 
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Fig. 5. Comparison of Eq. [T] {dashed) and our results {solid, 
dotted) for the Bonnor-Ebert sphere model. 



x[ ^/^(Qr ■ aiM(4;r7,/c)]2(- 



1 + 



"^drag,g 
"^drag.em 



-f. 



(7) 



3. Polarization 

3.1. Rayleigh polarization reduction factor 

The Rayleigh polarization reduction factor R is a measure of 
the imperfect alignment of dust grains with respect to the mag- 
netic field (Greenberg |1968t see also Lee & Draine 1985 ). 
The degree of polarization is reduced when grains are not 
aligned with the magnetic field. In the case of radiative torques, 
smaller grains are not aligned but larger grains are. The mini- 
mum aligned grain size, Oaig, is derived from Eq. |7] such that 

fleff = flalg when Wrad/Wr > 5. 



R = 



The polarization reduction factor is 

X"™' Cr3nn{a)da 
X""" Ci-ann(fl) da 



(8) 



where Cran is the polarization cross section of a randomly 
aligned grain, n{a) the grain number density, a the grain size, 
flmin the minimum size of the grains, and Umax the maximum 
size. This integral is, of course, sensitive to the limits and the 
form of the grain size distribution (see Fig.[TJ. In the following, 
we use flmin = 0.005 yum and the Li & Draine (1200 II I grain size 
distribution. 

3.2. Polarized thermal dust emission 

The polarized thermal dust emission is calculated following 
the formalism in Fiege & Pudritz (2000). Self-absorption and 
scattering can be neglected at submillimeter wavelengths. The 
Stokes Q and U components are equal to the integrals 



al^ cos Itfr cos y ds 



ah sin lip cos y ds 



(9) 

(10) 



where c is a coefficient of the particle properties to be defined 
later, the angle between the projection of B on the plane of 
the sky and the north, and y the angle between the local B vec- 
tor and the plane of the sky. The dust emission intensity at the 
given wavelength, is obtained from radiative transfer calcu- 
lations. 

The polarization angle is given by 



tan 2x ■ 



q 



and the degree of polarization P is 



P = 



with 



h ds , 



and 



Q'/,i(cos y - —) ds , 



(11) 



(12) 



(13) 



(14) 



where 2 is the intensity of the dust emission along the line of 
sight, and 1,2 is a quantity related to that, which is determined 
by the extinction cross section (see Eq.lO in Fiege & Pudritz 
l2500l) . 

The coefficient a is defined as 



a^RF 



pol 



(15) 



where R is the Rayleigh polarization reduction factor due to 
imperfect grain alignment, F is the polarization factor due to 
the turbulent component of the magnetic field, Cpoi is the grain 
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Fig. 6. Simulated polarization maps at 353 GHz with R calculated using Eq. [T](left frame) and Eq. |2l Polarization vectors are 
drawn for every third pixel and the scaling of the vector lengths is the same in both frames, representing the polarization degree. 
In the left frame the maximum polarization degree is 8.5%, while in the right frame it is 7.3%. The background image shows the 
logarithm of the total intensity at this frequency. 



polarization cross section, and Cian is the average cross section 
of a randomly aligned grain. The polarization cross section of 
an oblate grain is the difference between the main axis cross 
section and the minor axis cross sections of the grain (for more 
details, see Lee & Draine [T985l l. 

In our study, F - 1 because the three-dimensional mag- 
netic field is given by the MHD models, and we assume that 
the small-scale structure is resolved in the numerical solution. 
The ratio Cpoi/Cran needs to be fixed because we want to study 
the efifects caused by variations in R. To avoid an unreasonably 
high polarization degree, we choose Cpoi/Cian = 0.15. This cor- 
responds to the axial ratio of roughly 1.1 (see Fig. 1 in Padoan 
et al. l200TT i. 

4. Radiative transfer calculations 

To calculate the amount of polarized dust emission, we needed 
to carry out radiative transfer modeling, to (1) calculate the in- 
tensity and angular distribution of the radiation field, (2) deter- 
mine the resulting total intensity of dust emission, and (3) de- 
termine the dust temperature. Apart from determining the total 
emission, dust temperatures affect thermal rotation rates and, 
thereby, the minimum size of aligned grains. The intensity of 
the incoming radiation and its anisotropy affect the magnitude 
of the rotational torques, as shown in Eq. |7] 

The radiative transfer calculations are made with a Monte 
Carlo program (Juvela & Padoan I2003I I. where the contribu- 
tion of transiently heated, small grains is also solved. To de- 
termine the anisotropy of the radiation field within the model 
clouds, the photons entering each computational cell were reg- 
istered. The angular distribution was tracked by binning the 
photons into 12 bins that corresponded to angular discretiza- 
tion according to the Healpix scheme (Gorski et al. 120051 1. The 
anisotropy of the radiation field was estimated by calculating 
the vector sum of the intensity over different directions. We 
made the simplifying assumption that the radiative torque efli- 



ciency is direction invariant and that opposing directions can- 
cel out. The remaining part was our anisotropic radiation com- 
ponent. When the anisotropic component was subtracted from 
the total radiation, we were left with the isotropic component. 
Based on Table 4 in Draine & Weingartner (1 996) , we assumed 
that 2p° = 0.1 X Qr, where Qr is the radiative torque efficiency 
for anisotropic radiation, shown in Fig. |2] We used the dust 
model discussed in Li & Draine ( I2001I I. 

Finally, we calculated an example where we took the di- 
rection dependence of the radiative torque efficiency from Fig. 
17 in Hoang & Lazarian (i2008bj . We normalized the unidi- 
rectional curve to 1 when the radiation vector and the mag- 
netic field vector are parallel, and used the resulting curve as 
our direction-dependent efficiency factor. The efficiency fac- 
tor drops by more than a factor of ten when moving towards a 
perpendicular orientation of the vectors. We considered neither 
multipoles nor isotropy in this simplified example, but multi- 
plied each beam direction by the efficiency factor, and calcu- 
lated the vector sum of our anisotropy. This anisotropic radia- 
tion intensity was then used to calculate the aligned grain size. 

5. The model clouds 

We examine briefly one spherically symmetric cloud model for 
which the density distribution is obtained from the Bonnor- 
Ebert solution of hydrostatic equilibrium (Bonnor 1956 Ebert 
I1955I I. The mass of the cloud is 3.7 Mq and, for a stability pa- 
rameter value of ^ = 7.0, the visual extinction to the cloud 
center is Ay ~15'". When the total opacity of the cloud is fixed 
the exact value of the parameter ^ is not important. For a given 
distance from the cloud surface, measured in Ay, the local in- 
tensity and anisotropy of the radiation field depend only weakly 
on the true shape of the density profile. 

Our main interest lies in the study of inhomogeneous, 
three-dimensional clouds. We revisit the model C discussed by 
Pelkonen et al. (I2007I I. It is based on the results of numerical 
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Fig. 7. The relation between polarization degree and total intensity for selected cores of the polarization maps in Fig.|6] Cyan 
dots show the old results (Eq.lTJ, assuming constant y, while black dots show the new ones (Eq.|7]l with spatially varying y. 



simulations of highly supersonic magnetohydrodynamic turbu- 
lence, completed for a 128^ computational mesh with periodic 
boundary conditions (for details, see Pelkonen et al. I2007I I. 

In addition, we study three cloud cores that are taken from a 
new simulation of supersonic and super- Alfvenic MHD simu- 
lation. The MHD simulation was run on a mesh of 1000^ com- 
putational cells with the Stagger Code (Padoan et al. [2007 ), 
re-meshed to 1024^. In the used snapshot the rms sonic Mach 
number is 8.9 and the rms Alfvenic Mach number is 2.8, so 
that the turbulence is still super- Alfvenic even with respect to 
the rms Alfven velocity. From the 1024-^ volume, we selected 
three dense regions ('clumps'), each 128^ cells in size (see 
Fig. |3]l. The average density of the 1024^ cloud was scaled to 
n(H2)=150cm"-' and its linear size was set to be 6pc. This cor- 
responds to an average visual extinction of ~3™. In the selected 
core regions the visual extinction is higher, reaching peak val- 
ues of 10.9, 23.7, and 21.8 magnitudes in the three cores, re- 
spectively. The radiative transfer is solved locally for the sub- 
grids of 128"^ grid points rather than the full 1024^ mesh. 

6. Results 

6.1. Spherically symmetric cloud 

We start by looking at grain alignment in the case of the spher- 
ically symmetric cloud described in Sect. |5] Figure |4] shows 



the anisotropy factors at different wavelengths. The intensity 
of the short wavelength radiation drops sharply with Ay, but 
the anisotropy increases since the remaining intensity comes 
preferentially from smaller and smaller space angles. At longer 
wavelengths, the dust emission starts to dominate over ISRP 
and the direction of the anisotropy is towards the surface, al- 
though the effect on effective y is small. At the center of the 
spherically symmetric cloud the anisotropy is by definition 
zero. 

Figure|5]shows the comparison of Eq. [Hand our results with 
constant y - 0.7 and y calculated from the radiative transfer 
calculations. The number of grains larger than 0.58 fim in Li 
& Draine (12001!) dust model is negligible. Thus, Fig.[5]shows 
that according to this calculation, there is essentially no align- 
ment above 3 magnitudes, while in the Cho & Lazarian model 
the limit is 5 magnitudes. Larger grains, if they exist, could be 
aligned by radiative torques at higher Ay. In an inhomogeneous 
cloud there might also be some sightlines of lower Ay, which 
produce more efficient radiative torques. 

6.2. Three-dimensional MHD models 

The dust polarization was studied in three-dimensional models 
described in Sect. [5] Figure[6]shows the calculated polarization 
maps at 353 GHz for our old 3D model C, used in Pelkonen 
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Fig. 8. Simulated polarization maps of the high resolution model at 353 GHz. Polarization vectors are drawn for every fifth pixel, 
the length corresponding to the polarized intensity. The background image shows the logarithm of the total intensity. The scaling 
of the background and the length of the vectors are the same in all frames, the longest vector being 0.11 MJy/sr. 



et al. (120071) . The polarization degree has been reduced every- 
where, as expected from Fig. |5] but particularly in high den- 
sity regions. This reduction in the polarization degree is more 
readily apparent in Fig.|2l where the Pjl slopes of most of the 
selected cores are steeper with the new calculations. This is be- 
cause in a clumpy medium, starlight can penetrate the cloud 
from many different directions rather than in a clearly defined 
'brightest' direction. Thus, the anisotropy will be averaged out 
somewhat, leading to weaker radiative torques and lower po- 



larized emission. In the cores themselves the anisotropy is even 
less, so the difference with an invariant j - 0.7 is greater, re- 
sulting in a steeper slope. 

Figures [8] and |9] show similar plots for our new high reso- 
lution model at 353 GHz. In Fig. [8] the vectors are drawn for 
polarized intensity, not for polarization degree, which is why 
long vectors are superimposed on the bright, high density re- 
gions. In Fig. |9l Core 1, when viewed from the z direction, 
does not seem to show a clear decrease in the polarization de- 
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Fig. 9. The relation between polarization degree and total intensity for the polarization maps in Fig. [8] In order to avoid cluttering 
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Fig. 10. The anisotropy factor at A - 0.55 nm in an xy 
slice. The slice is the mid-plane of the high density model 
for Core 2. The contours mark the density at log«//, = 
3.5,4.0,4.5,5.0,5.5,6.0. 



gree, unlike the other examples. This is because Core 1 is less 
dense, and the maximum visual extinction through the core in 
the z direction is only 3.9 magnitudes. This also explains why 
the other two relations for Core 1 are not as steep as for Cores 
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Fig. 11. The relation between polarization degree and total in- 
tensity for Core 2 (see Fig.|9|l where the density has been multi- 
plied by ten. The figure shows the polarization degree observed 
towards the z direction. The polarized dust emission is calcu- 
lated for doubled (approximate calculation, see text) and nor- 
mal grain sizes, and for a mixed distribution where larger grains 
are found only in dense regions (see text). 



2 and 3, because radiation can penetrate the core much more 
easily. 
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Fig. 12. The gas density, (solid line) and the polarized intensity (dashed line, at 353 GHz) along one line of sight through 
each core. The plots represents a line of sight through the intensity maximum of each map in Fig. |8] For Core 2, the polarized 
intensity is also calculated for doubled grain sizes (dotted line). The number in the lower right corner of each plot is the visual 
extinction along this line of sight. 



In dense clumps, grain coagulation can cause the grains 
to grow, shifting the size distribution to larger grain sizes. We 
took Core 2 and multiplied the density by ten to obtain a dense 
cloud. Figure [TO] shows the anisotropy in a slice in the mid- 
plane, as well as density contours. In Fig.[TTl we show the po- 
larization degree for three different dust size distributions for 
this high density model. We first calculated the grain alignment 
and dust emission using normal Li & Draine (I2001I I dust. The 
polarization degree is rather low and declines over the com- 
plete intensity interval. Secondly, we calculated the alignment 
by assuming that all our dust grains were double their original 
size. The shape of the radiative field inside the cloud should 
be closely the same, and we were only interested in the po- 
larization degree and scaled intensity rather than absolute val- 
ues of intensity. Thus, we used the same Oaig as for normal 
dust and simply used the doubled grain size distribution shown 
by Eq. |8] resulting in more efficient polarization. There is a 
clear plateau in low intensity, before the grain size needed for 
alignment becomes larger and the polarization degree begins to 
drop. Finally, we investigated a model in which large grains ap- 
pear only in dense regions. The two dust populations, one with 
normal grain size distribution and another with twice as large 



grain sizes, were calculated from the formulas: 

= l+e-log«,/0.3' (16) 

X„ = l-Xd, (17) 

where Xd and X„ are the abundances of double sized and nor- 
mal dust, and ns = «H/10^cm"-'. The result, as expected, is 
between the two previous cases. There is a plateau, but it ends 
sooner than in the large grains case. However, towards the high 
intensity, the slope is becoming less steep, due to the grain 
growth in the densest clumps. 

Where does the observed polarized emission originate? The 
thermal dust emission tends to trace the density structure. High 
density means more emitting dust and thus higher emission, 
although the temperature is anticorrelated with density and thus 
complicates the picture. However, Figs.[T]and |5]hint that above 
a couple of Ay only a fraction of the grains are aligned, and 
thus the polarization degree drops sharply. Figure[T2]shows the 
polarized intensity and the gas density in each cell along lines 
of sight through each of the cores as seen from three different 
viewing directions. The effect of shadowing by even modest 
density structures is striking, as is the absence or reduction in 
the polarized emission from the densest structures. In Core 2, 
the polarized emission is also plotted for the doubled grain size. 
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Fig. 14. The RMS difference between the observed polarization angles and the mass-averaged (dashed line) and volume-averaged 
(dotted line) magnetic field directions, for all cores and viewing directions. The observed polarization angles are calculated with 
normal Li & Draine (2001) grain size distribution. For Core 2 polarization was also calculated with doubled grain sizes, and 
compared with mass-averaged ((thick solid line) and volume-averaged (thick dash-dot line) magnetic field direction. 



In this case, the polarized emission originates deeper inside the 
core. However, because of the high local density, the density 
maximums in the y and x directions still do not have aligned 
grains. The z direction exhibits significant polarized emission 
from the density maximum even though the visual extinction 
along that line of sight is greater. In Core 1, the magnetic field 
is oriented almost towards the observer in the z direction, which 
is why the polarization is weakened. 

Figure [13] shows the histogram of both the intensity and 
the fraction of polarized intensity in different density bins for 
Core 2 seen along the z direction. Due to the summing of the 
contribution of the individual cells in each bin, we need to com- 
pare the total intensity contribution with the polarized emission 
contribution. The high density cells contribute less and less to 
the polarized emission, while the contribution of the small den- 
sity cells is enhanced. With doubled grain sizes, the more ef- 
ficient grain alignment is clearly evident because the polarized 
emission still comes from higher density regions, before finally 
falling, too. 

Do polarization vectors trace the magnetic field with any 
degree of certainty? Above a threshold, which depends on the 
local density as well as the geometry of the cloud, the radiative 



torque alignment is no longer capable of organizing the grain 
orientation. Thus the emission from shielded region is unpo- 
larized, giving the observer no information about the magnetic 
field within it. 

Figure [14] shows the rms errors between the observed po- 
larization angle, Xobs, and the volume- and the mass-averaged 
magnetic field angles, xb.yoI and XB,mass- For low Ay, XB,mass 
and Xobs agree well, because there are no dense clumps along 
the line of sight where the radiative torque alignment would 
be inefficient. Temperature differences between warm diffuse 
regions and cold denser regions do introduce some discrep- 
ancies caused by the different weighting in emission and in 
density. For higher Ay, the situation changes; ;^B_inass traces the 
magnetic field directions in the dense clumps, while the influ- 
ence of dense clumps on Xobs is reduced by the weakened ra- 
diative torque alignment. This is because the radiative torque 
alignment continues to be effective on the larger grains, and 
the dense clumps continue to contribute to the polarized emis- 
sion. Any emission or density variations are not taken into ac- 
count by XB,voh and thus the difference between observed and 
volume-averaged angles depends greatly on the density profile 
and the local magnetic field orientation. If the density profile 
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Fig. 13. A histogram showing the contribution of different den- 
sity cells to the total intensity (solid line), and to the polarized 
intensity for normal grain sizes (dashed line) and doubled grain 
sizes (dotted line). This is for Core 2 and the local polarized 
emission as seen from the z direction. The total and polarized 
intensity are shown as sums of all emission from the cells in 
that bin; thus the few high density cells have a low contribution 
to total intensity. All intensities are at 353 GHz. 



is flat or the magnetic field is ordered apart from the densest 
clumps, thenxB.mi ~ Xobs- 

At the highest Ay, there are only a few independent struc- 
tures (see Fig fTSl i. and the different magnetic field geome- 
tries play a crucial role. Taking Core 2 as an example, there 
is one structure along the z-axis. In that structure, two clumps 
along the line of sight have their magnetic fields aligned anti- 
parallel. Their combined signal is mostly unpolarized, and thus 
Xobs ~ XB,m^ss ~ Xb.yoI- That is why the rms remains flat even 
to high Ay. The same would happen if both clumps had a mag- 
netic field direction close to;t'B,voi, in which case the contribu- 
tion of the clumps - whether total (^B,mass) or partial (xobs) - 
would only reinforce the orientation from the diffuse regions. 
Along the y-axis, the density profile is quite varied, as is the 
magnetic field orientation, and thus none of the angles agree 
at high Ay. Finally, along the x-axis, direction of the mass- 
averaged magnetic field is determined mostly by one high den- 
sity clump, where the magnetic field direction spins around. 
Thus, Xoh>. + ;rB,mass but;t'obs ~ XB,^a\- As always with polariza- 
tion, one cannot ignore the geometrical effects. 

From Fig. [14] it can be seen that an rms difference of 
10-30 degrees is common for sightlines of higher extinction 
than a few Ay. This difference arises because the grains are 
not aligned in the dense regions of the cloud, and thus the 
dense regions do not contribute to the observed polarization. 
This may cause problems for the Chandrasekhar-Fermi method 
(Chandrasekhar & Fermi |1953l l and in attempts to infer core 
formation processes based on the relative alignment of B and 
core geometry. 

Despite the above caveat, the Xobi and XB,m3.%'i do seem to 
agree until the extinction exceeds a couple of Ay. If grain size 
increases (Core 2, solid line), Xob% and XB,msssi can remain in 
a good agreement even to Ay ~ 10™, as predicted by Cho & 




Fig. 15. The gas density, n^^, (solid line) and the polarized 
intensity (dashed line) along one line of sight through Core 2. 
The polarized intensity is also calculated for doubled grain size 
(dotted line). The number in the lower right corner is the visual 
extinction. 



Lazarian ( I20051 l. This is because the radiative torque alignment 
continues to be effective on the larger grains, and the dense 
clumps continue to dominate the polarized emission. Figure [TSl 
shows an example along one line of sight where the doubled 
grain size allows the dense clump to be probed. 

Finally, we examine the effect of the dependence of the 
radiative torque efficiency on the angle between the radia- 
tion field and the magnetic field. We limit the study to the 
anisotropic component of the radiation field. In Fig. [161 we can 
see the clear need for larger grains because of the weakening of 
the radiative torques at larger angles. Some isolated cells show 
the opposite effect, where the summed radiative torque is ac- 
tually stronger due to less efficient opposing torques. Usually, 
there is a drop of a factor of two in the polarization degree of 
the denser regions. 

The weakening of the radiative torques means that we 
miss more of the polarized emission from dense clumps in the 
clouds. This limits our ability to trace the mass-averaged mag- 
netic field direction along sightlines of high Ay. In the normal 
grain size case, our ability to trace the magnetic field dropped 
from Ay ~ 3™ to 1™. However, this is a complex issue. In 
Fig. [141 we can see examples where magnetic field is trace- 
able over Ay - 10" even for normal grains, because of the 
orientation of the magnetic field. Also, the Ay limit depends 
on the direction the core is viewed from. For example. Core 

2 gives quite different results for each different viewing di- 
rection. Changing our alignment parameter, Wmd/'^r = 5, to 

3 (see Discussion) would also improve the alignment, helping 
to counteract some of the weakening of the radiative torques. 
Grain growth is of critical importance. Based on our current 
results, we conclude that Ay = 10" is likely to be an upper 
limit, and polarized emission provides reliable data on mag- 
netic fields only up to a few magnitudes in Ay. 

7. Discussion 

One critical parameter for grain alignment by radiative torques 
is the grain size. Figure [T4| shows that in the case of a diffuse- 
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Fig. 16. Upper panel: The aligned grain sizes for direction- 
dependent radiative torques, flaig.dir, and direction-invariant ra- 
diative torques, aaig- This is for Core 2, taking into account 
only the anisotropic component. The solid line shows where 
flaig.dii = «aig- The gaps in the grain sizes is due to the size grid 
used in the calculations. Lower panel: The relation between po- 
larization degree and total intensity calculated using the above 
grain sizes and seen along the z-axis. 



medium dust model the observed polarization vector no longer 
traces the mass-weighted magnetic field direction in clumps of 
extinction higher than a few Ay. However, if the grain size dou- 
bles, the polarization vectors can carry information about far 
denser clumps. This is highly dependent on the geometry of 
the cloud, as can be seen in Core 2, where, along the y-axis, we 
reach only Ay ~ 7™, while along the x-axis we reach ~ 10™. 
A doubling of the grain size is not improbable: the sub-mm 
emissivity of dust has been observed to increase significantly 
in clouds that have a visual extinction of only a few magni- 
tudes (Bernard et al. |T999l Stepnik et al. |2003l l. This has been 
interpreted as evidence of grain growth and also of the appear- 
ance of dust aggregates and fluffy grains of low packing den- 
sity. Whittet ( 2007 ) found evidence in a globule of dust growth 
to mean sizes 50 % - 100 % larger than in the diffuse ISM, even 
along a sightline of Ay - 1.1. 

For radiative torques the important question is the growth 
of the largest grains. Ossenkopf & Henning C1994) presented 
several models of grain growth with thin and thick ice man- 



tles, and coagulation. They were able to explain the observed 
changes in dust opacity. In their model of thick ice mantles, 
applicable to cold dense cores, the ice volume fraction was a 
factor of 4.5 above the original dust volume, causing a signif- 
icant change in the size distribution. But since the ice mantle 
thickness is almost constant irrespective of the grain size, most 
of this increase in size happens in the smaller grains. Grain 
coagulation is needed to make the large grains larger, and this 
may take longer than the free fall time of a cloud (Ossenkopf & 
Henning [T994l Chakrabarti & McKee 2005). Hence, the mag- 
netic field within the youngest cores may not be accessible, 
while it could be probed inside older cores and filaments that 
survived significantly longer than their free-fall time. 

Many observations of the dust polarization at sub-mm 
wavelengths focus on clouds that are rather luminous in dust 
emission (e.g. Davis et al. 120001 Matthews & Wilson 120001 
120021 Henning et al. 120011 Wolf et al. l2003l Lai et al. l2003l 
Crutcher |2004| l, many with internal sources. Thus, the compari- 
son of our results derived from a MHD modeling with cold dust 
and without internal sources is not straightforward. However, it 
is interesting to note that if the intensities are normalized by the 
maximum intensity of a given core, the modeled Pjl -relations 
show qualitatively similar behavior to the observed ones, drop- 
ping from 10% to 1% in polarization degree. 

Cho & Lazarian ( I2005I I chose their alignment parameter, 
Wiad/wr - 5, and derived Eq.[T]from that assumption. Lazarian 
& Hoang ('2007 ) presented a simple toy model of a helical grain 
to study the basic properties of radiative torque alignment, and 
in their later work (Hoang & Lazarian 2008a) found that when 
Wrad/wj- = 3, a significant fraction of the grains are aligned. 
Bethell et al. ([2007) decided to use w,ad/wr = 3 in their own 
study of the anisotropy of the radiation field inside a clumpy 
cloud. We chose to use the higher value to be able to compare 
the results to our previous paper (Pelkonen et al. 120071) . The 
lower alignment parameter leads to a more effective alignment 
of grains. This does not dramatically change the results of this 
study, since, as seen in Fig |2l Eq. |7] is very sensitive to the 
grain size. As a result, we need only slightly larger grains than 
BetheU et al. dlOOTl l. 

Direct comparison with the results of Bethell et al. ( I2007I I 
is further complicated by the different cloud models and 
grain size distributions. However, with their mass-averaged 
anisotropy factor y - 0.34, we would not expect isotropic ra- 
diative torques to play any significant role. In our high-density 
model that showed similarly high anisotropy, isotropic contri- 
bution was negligible. In the low-density case, anisotropy was 
low and thus isotropic contribution was noticeable as a small 
enhancement of polarization degree. Despite the setup differ- 
ences with the study of Bethell et al. ( I2007I I. we have obtained 
qualitatively similar results and hence agreement with their 
conclusions. However, we would like to add one note of cau- 
tion to their statement about invariant anisotropy having only a 
moderate effect on the polarization degree. When looking at the 
emergent polarized emission along a line of sight, it is obvious 
that often the very densest clumps do not contribute at all. This 
is due partly to the weakening of the radiation field, but also 
to the smaller anisotropy inside the core. Adopting an invariant 
anisotropy would decrease the contribution of the shell around 
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the core, where the anisotropy is highest, and might cause the 
grains in the dense core to be aligned. Hence, while the big pic- 
ture of the Pj I relation might look similar, the magnetic field 
that is sampled in the model might be quite different. 

8. Conclusions 

We have examined the anisotropy of radiation in spherically 
symmetric cloud models and in inhomogeneous 3D cloud mod- 
els. In our ID cloud model, the anisotropy was well below 
0.7 (Fig. resulting in a weaker radiative torque alignment 
(Fig.EJ. 

In a 3D cloud, the line of sight Ay can be high, while the 
effective Ay can be low due to inhomogeneity. These low Ay 
sightlines can reduce the anisotropy because of averaging and 
thus reduce the polarized signal, as also noted previously by 
Bethell et al. (120071) . Total anisotropy is sensitive to the ge- 
ometry of the cloud, but also to the Ay of the cloud. The P// 
relations show a qualitatively similar behavior to the observed 
ones. 

The 3D model and already the ID case show that radia- 
tive torque alignment does not operate well inside cloud cores, 
where Ay is larger than a few magnitudes. Thus, we should be 
careful not to interpret the polarization that we see from a dense 
core as originating in the core, but rather from the more diffuse 
regions on the line of sight, where the ISRF is able to align 
the grains via radiative torques (Fig.fTZt. These results suggest 
that one must be careful with models of core formation. We in- 
tend to study the impact on the Chandrasekhar-Fermi method 
in future papers. 

The growth of grains in the clumps could result in a more 
efficient alignment of grains even in the denser clumps, allow- 
ing them to continue contributing to the observed polarized 
emission. Thus, if grain growth occurs in the clumps, the ob- 
served polarization vectors may trace the magnetic field lines 
possibly up to Ay -10 magnitudes, in agreement with the pre- 
diction in Cho & Lazarian (I2005I I and the results of Bethell et 
al. (|2007| |. However, prestellar cores are often found to evolve 
on dynamical timescales, suggesting that larger grains would 
not have time to coagulate, and polarization vectors would not 
trace the magnetic field within the cores. Furthermore, using 
radiative torques that account for the angle between the radia- 
tion direction and the magnetic field, the polarization is again 
reduced. As a result, dense cores would need even larger grains 
for their inner magnetic fields to be detectable. Even with grain 
sizes at twice the size, we are likely to trace magnetic fields 
only up to a few magnitudes in Ay. 
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